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ABSTRACT 

We investigate the formation of star clusters in an unbound GMC, where the supporting 
kinetic energy is twice as large as the cloud’s self-gravity. This cloud manages to form a 
series of star clusters and disperse, all within roughly 2 crossing times (10 Myr), supporting 
recent claims that star formation is a rapid process. Simple assumptions about the nature of the 
star formation occurring in the clusters allows us to place an estimate for the star formation 
efficiency at about 5 to 10%, consistent with observations. We also propose that unbound 
clouds can act as a mechanism for forming OB associations. The clusters that form in the 
cloud behave as OB subgroups. These clusters are naturally expanding from one another due 
to unbound nature of the flows that create them. The properties of the cloud we present here 
are are consistent with those of classic OB associations. 
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1 INTRODUCTION 

Practically all star formation is thought to occur i n clusters that 
are e mbedded in giant molecular clouds (GMCs) iUada & Ladal 
l2003h . This suggests that a clustering environm ent, where multi¬ 
ple objects co mpete for a common gas reservoir iZinneckeJ 1 9821 : 
|LMS^|T99a), plays an important role in early stages of protostel- 
lar evolution, s uch as dictating the form of the stellar in itial mass 
functi on (IMF) jBo nnell et a lll997l : lBonnell et all200ll) . Further¬ 
more, lElmegreenI i2Q0(J) has collected observational evidence sug¬ 
gesting that star formation is a rapid process, occurring on roughly 
the crossing time of the region at a variety of scales. Not only does 
he propose that the star formation in a typical GMC occurs within 
~ 4 Myr (approximately the crossing time for standard GMC) but 
that the cloud’s dispersal occurs within a few crossing times, or < 
10 Myr. 

The combined implication of these observations is that 
star formation occurs quickly and in groups and that the sites 
of star formation disperse quickly. Our proposal in this pa¬ 
per is that this is possible if GMCs are dynamically un¬ 
bound objects, with the internal turbulent energy greater than 
that of the cloud’s self-gravity. This follow s from the work of 
IVazQuez-Semadeni. Passot & Pouqu^ ^1995^ who showed that 
transient (unbound) GMC sized objects can be formed from flows 
in the ISM. We also find that unbound GMCs may provide a nat- 


* E-mail: pcc@st-and.ac.uk 


ural mechani sm for the creation of OB associations, a notion first 
suggested bv lAmbartsumianI il9.58l) . 

In the rest of this first section we discuss the ideas behind rapid 
star formation and the dynamical state of GMCs. We also include 
in this section a discussion of OB associations. In section|^we de¬ 
scribe the details of the simulation and section0follows the general 
evolution of the GMC. In section|3we give estimates of the star for¬ 
mation efficiency in the GMC based on some simple assumptions. 
In section 13 we highlight the similarities between the simulation 
and the general structure in an OB associations. A summary of the 
paper’s main conclusions can be found in section|3 


1.1 GMC Lifetimes and Rapid Star Formation 

Until the last decade or so GMCs were generally believed to 
be long-lived structures, w ith some estimates of ages reach¬ 
ing as high as 10^ Myr Solornon. Sander s & Scovillj 1 19791 : 
IScoville. Solomon & Sand^ 1l97^ ~ lscoville & Hershi Il979l) . It 
was generally believed that the chemistry of turning atomic species 
into molecules would require millions of years befo re an objec t 
like a GMC would be detectable via its CO abundance iJurJl97^ . 
One also had the problem that the CO mass in the galaxy, cou¬ 
pled with estimates of the star formation rate, suggested that GMCs 
had to live for tens of millions of years if the star formation ef- 
ficiency was to remain at the observed level of a fe w percent 
iZuckerman & Evans! 197^IZuckerman & PalmeJl974h . 

Recent observations of embedded clusters tend to suggest that 
the whole process of star formation, including GMC formation 
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and dispersal, occ urs on roughly the crossing time for the region 
iElmegre^bOOCh . Not only do most molecular clouds in the so¬ 
lar neighbourhood contain signs of star formation in the form of 
clusters, but the age determination of thes e clusters suggest s they 
are very young, typically less than lOMyr jHartman t]l200(t . This 
suggests that star formation occurs quickly in GMCs after their for¬ 
mation. The fact that clusters with ages greater than ~ 5 Myr are 
seldom associated with molecular gas, sugge sts that clouds disperse 
quickly iLeisawitz, Bash & Thaddeusll989h . 

In the original cloud lifetime proposition, it was assumed that 
all of the CO observed in the galaxy was associated with molec¬ 
ular hydrogen involved in star formation. We now realise that the 
vast majority of the gas that comprises a GMC is never involved 
in the star formation process. In fact the star formation efficiency 
in GMCs is only a few percent. The reason behind this lies with 
the fact that little of the cloud is actually dense and bound enough 
to turn into stars in th e cloud’s lifetime iPadoan|[l995tlitotmannl 
ll998tlZinneckeJ2002h . Also, if GMCs are short-lived features then 
there is little time for the more tenuous parts of the cloud to get in¬ 
volved in the star formation via accretion. 

The old GMC model also required the cloud to be supported 
and in virial equilibrium, since that would permit them to remain as 
coherent structures for as long as was necessary. This support pres¬ 
sure had to be in the form of non-thermal kinetic energy, such as 
turbulence iLarsonlll98lh . since the thermal energy component of 
these clouds is typically very small. To counteract the gravity on the 
large scales however requires motions which are supersonic, and it 
is known that these quickly damp in shocks iMac Low et all 19^ 
IStone. Ostriker & Gammiell9981) . even in the presence of magnetic 
fields. Thus the bound GMC model requires some method of con¬ 
tinually driving the turbulence on the large scale. These driving 
mechanisms are not necessary in the short cloud lifetime model, 
and there also is no need to assume that the clouds are in virial 
equilibrium. GMCs can therefore exist in a variety of dynamical 
states^ _ 

iHever. Carpenter & Snelll i200lh have examined the stability 
of molecular clouds in the outer galaxy and come to the conclu¬ 
sion that most clouds are indeed globally unbound by their internal 
motions. They also point out the difficulty in producing mass es¬ 
timates (which a great number of papers on the subject of GMCs 
pass over) and note that even mass estimates determined via CO 
measurements (both ^^CO and ^®CO) assumes at some stage the 
cloud is bound. Although they do find the clouds approach dynam¬ 
ical stability at large masses (> IO^Mq), there is still considerable 
scatter in the data, suggesting that there is no typical dynamic state 
for GMCs. _ _ 

IPringle. Allen & Lubowl i200lh have shown that it might also 
be possible to build GMCs by accumulating very low density hy¬ 
drogen gas, already in a molecu lar state. Their stu dy came in re¬ 
sponse to the ideas presented bv lElmegreenI i200dl . in an attempt 
to provide a new mechanism for GMC formation that can occur 
quickly. They point out that it is quite possible that a large fraction 
of the interstellar medium may be in molecular form, but either sim¬ 
ply too low a density to be detectable by current methods or too far 
away from illuminating sources. The GMCs are then formed from 
large scale shocks, from spiral arm passage or feedback from high 
mass stars, such as winds and su pernovae. This cloud formation can 
occur within a few million vears. lpringle et ^ J200il) also point out 
that GMCs are probably not in virial equilibrium, and note that their 
wind-swept appearance suggests that they are anything but. 

The simulation that we present here draws on the above stud¬ 
ies for motivation. We assume that large scale flows are able to 


create an unbound GMC in a few mi llions years. Instea d of be¬ 
ing contained by external forces (e.g. iHever et al]l200lll . we as¬ 
sume that the cloud is free to expand into the ISM. Thus the 
flows that created the cloud are assumed to have been used up 
in its formation. Since the cloud is assumed to be short lived 
and not quasi-static, there is no need for the internal turbulent 
energ y, which will dissipate on the crossing time, to be replen¬ 

ished jBa ll esteros-Paredes. Hartmann & VazQuez-Semadeniil999l: 

lElmegr^nl2'oo'cr 

1.2 The Origin of OB Associations 

OB associations are historically identified simply as extended 

i roups of OB stars , having diameters of tens of parsecs 
Ambartsumian|[l95^ . Furthermore they are rather more diffuse 
than open c lusters, with th e mass density of OB t y pe stars at ~ 0.1 
Mn pc~ ^ iBlaauwl ll964l: lAmbartsumianI 1 195.51 : lGarman^ll994l: 
iLada & Lad^l200^ . It was found that these associations contain 
considerable s ubstructure which are referred to as ‘OB subgoups’ 
iBlaauwll964l . These subgroups are unbound from one another as 
was deduced f rom their expansion about the centre of the region 
^Blaau^^^l9■52^1 . Some regions or ‘subgroups’ are shown to be asso¬ 
ciated with molecular gas. In general these regions are not coeyal 
but can exhibit a sp read of ages be tween the subgroup population 
as large as 10 Myr iBlaauwlll964h . The fact that OB associations 
are yery young, with some of the subgroups possessing ages of the 
order of a millions years, suggests that unbound nature of the sub¬ 
groups from one another is primordial. 

The relationship between OB associations and other types 
of clusters found in the galactic disc, such as open clusters 
and embedded clusters, is still rather unclear. The OB associa- 
tions do howeyer haye a classic theory regarding their formation. 
lElmegreen & Lad^ il977h proposed that OB associations form via 
triggering, prompted by the ionised regions produced by previous 
generations of OB stars. In this manner, the star formation is self 
propagating, with one generations of OB stars triggering the for¬ 
mation of the next. Since the shocked layer in which the new group 
of OB stars forms is moving away from the older OB stars, at a 
few kms“^, the new group is unbound from its parent group. The 
region then naturally has the dynamics of the observed OB groups. 
Motivation came from observations of stars forming at the bound- 
aries of molecular clouds and HII, such as NGC75 38, M17 andM8 
iHabing. Israel & de Jongll972l : ll^da et alll97fll . 

The issue is complicated however whe n one considers the de¬ 
tailed stellar population of OB associations iGarmanvl 1 994lBro^ 
1200 ih . In the self propagating model, OB type stars form in the 
shocked layers where conditions are naturally more suited to form¬ 
ing high mass stars. Low mass stars form spontaneously in the rest 
of the cloud. Thus the model assumes a two step formation process 
whereby low mass stars and high mass stars are formed by differ¬ 
ent mechanisms and in physically separated locations. The IMF of 
the OB associations however do not exhibit this feature and gener¬ 
ally possess the stan dard field star IMF, at least within the Salpete r 
range (e.g. Sco OB2, Ide Geusll992l : lfteibisch & ZinneckeJl999h . 
Since up to nearly 90% of star formation is thought to occur in 
embedded clusters, with a fi eld star IMF and pri mordial mass seg¬ 
regation (for a discussion see lLada & Ladal2003h . it may be that the 
formation of OB associations has more in common with standard 
clus tered star formation. 

_ [Bonnell. Bate & Vind i20Q3h and iBonnell. Vine & Batel 

i200^ have modelled cluster formation in a turbulently supported 
cloud. They modelled a IOOOMq molecular cloud that was 
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initially supported against collapse by a turbulent velocity field. 
It was found that the dissipation of the large scale supersonic 
flows produced a numher of distinct subclusters. Each subcluster 
contains at the core a massive star. The subclusters were mass 
segregated and each had a protostellar population consistent with 
that of the observed field star IMF, both of which are the result of 
competitive accretion. Since the cloud was initially bound, even 
more so after the dissipation of the turbulent energy, the whole 
system of subclusters are themselves bound to one another. They 
quickly merge within roughly 0.5 Myr (roughly twice the free-fall 
time for the original cloud). If this merging process was to occur 
on large scales, such as a whole GMC, one would never be able 
to form OB associations. The massive stars at the centres of the 
subclusters would find themselves in one large cluster. 

Our proposal in this paper is that OB associations are just a 
series of clusters that form in unbound GMCs. The expanding cloud 
produces a series of clusters that are unbound from one another due 
to the fact that the flows that form them are also unbound from one 
another. The clusters, which become OB subgroups, then simply 
expand away from their mutual centre of mass along with the gas 
from the cloud, rather than merge into a single cluster. Thus only 
one star formation mechanism is at play here: clustered formation. 
The OB association therefore will have the universally observed 
IMF. 


2 DETAILS OF THE GMC SIMULATION 

The fluid was modelled using the Lagrangian pa rticle method 
of smoothed particle hyd rodynamics, or SPH iLucvl 1 19771 : 
iGingold & Monaeharll977ll . The smoothing lengths are variable in 
both time and space, with the constraint that there must be roughly 
constant number of neighbours for each particle, which is chosen to 
be roughly 50 (with a fluctuation from 30 to 70 neighbours). We use 
the sta ndard artificial viscosity suggested bv iGingold & MonaghanI 
il983h with a = 1 and /3 = 2. Gravitational forces are calculated 
using dipole and q uadrupole moments obtained via a tree structure 
iBenz et alJI 199(1) . which is also used to construct particle neigh¬ 
bour lists. The code has been parallelised by Bate using OpenMP 
and the simulation presented here was performed on the UK Astro- 
physical Fluids Facility (UKAFF). 

Our simulation starts with a uniform density sphere of molec¬ 
ular hydrogen of radius 20pc with a mass of 1 x 10® Mq . 
The gas is isothermal and has a temperature of lOK. These num¬ 
bers (mass, size and temperature) are typical o f those reported for 
GMCs in the solar neighbourhood jBlitjl99lh . We model the gas 
with 500,000 SPH particles and are thus able to accurately fol¬ 
low th e formation of self-gravitating regions down to a mass of 
2 OM 0 iBate & Burkeril9^ : IWhitworthll99il . The free fall time 
associated with this cloud, the time taken for the unsupported gas 
to collapse under gravity to a central point, is roughly 4.7 Myr. The 
cloud has an initial Jeans mass of 30.4Mq . We do not include any 
feedback processes, such as stellar winds and jets or the effects of 
massive stars such as ionisation fronts and supernovae. 

To model the turbulence, we support the cloud with a Gaus¬ 
sian random velocity field with a power spectrum of P{k) oc 
which is consistent with a velocity field with a Larson-type rela¬ 
tion of a oc LP'^ where cr is the velocity dispersio n and L is the 
length scale of the region iMvers & Gammiell999tl . At the begin¬ 
ning of the calculation the ratio of gravitational to kinetic energy 
is 0.5 (Ekin = 2Egrav). We stress that the turbulent kinetic energy 
is able to decay freely in this simulation since we include no driv¬ 


ing mechanism. The timescale for the energy decay is the crossing 
time iMac Low et alll998l: ls tone et alll998h which for the initial 
velocity field is t^r = 4.2 Myr, slightly less than the free fall time. 

The SPH code includes the modification bv iBate et alJil99^ 
which replaces dense, self-gravitating, regions of the gas with point 
masses, or ‘sink particles’. These sinks allow the code to model 
the dynamical evolution of accreting protostars, without integration 
time steps becoming prohibitively small. We set the sink particles 
to form at a density of 1000 times the initial density, with a sub¬ 
sequent accretion radius of 0.17pc. When a particle finds itself at 
the centre of a dense, bound and collapsing region it is turned into 
a sink particle and its 50 to 100 neighbours are accreted onto it. 
With the resolution used for this simulation the sink particles start 
with a mass of at least I 5 M 0 before further accretion. Therefore 
we cannot think of these point mass objec ts as ‘protostars’, as was 
the case ir jBate. Bonnell & BrommI i2003l) . but instead assume that 
they represent ‘proto-clusters’. To prevent the ‘sink particles’ be¬ 
having as point masses in gravitational interactions, we smooth the 
sink-sink gravitational forces to a distance of Tmin = 0.2 pc in the 
form Fij = —Gmimj/(rij 4- rmin)^ between particles i and j. 

In the analysis that follows, we discuss the properties of star 
formation centres, or ‘SFCs’. These can either comprise of a single 
protocluster (or sink particle) or a coherent group of protoclusters. 
To identify SFCs where more than one sink particle is involved, we 
make use of the mass segregation that occurs naturally when the 
protoclusters intera ct in self gravitating groups (see for example 
iBonnell et a lJ2004h . First we sort all the protoclusters by mass. We 
then take the most massive protocluster and tag it and its fellow pro¬ 
tocluster neighbours within 0.5pc to be a members of ‘SFC 1’. We 
then go down the mass sorted list of protoclusters until we arrive 
at the next most massive protocluster that has not been associated 
with SFC 1. It becomes tagged as being a member of ‘SFC 2’. All 
the protoclusters within 0.5pc of this protocluster are now tagged 
as being members of SFC 2, unless they are already members of 
SFC 1. This process continues down the list of protoclusters un¬ 
til all have been assigned membership to a SFC. There is thus a 
resolution of 0.5pc which distinguishes one SFC from another. We 
find that the choice of 0.5 pc used in attributing memberships does 
not significantly affect the SFC population, since the protoclusters 
formed in the simulation are either well separated (and thus in iso¬ 
lation) or exist in dense groups. The radius of the SFC is given 
by the radius of the furthest protocluster from the centre of mass. If 
there is only one protocluster then the radius is simply the accretion 
radius, which is 0.17pc. 

One problem with trying to model turbulence in a numerical 
simulation of this type is that it is not always possible to resolve the 
velocity structure at all scales. Turbulence is assum ed to be hierar - 
chical, following a Larson-type relation of ct oc L°‘ iLarsonll98lh . 
In SPH, while a particle can have a kinetic energy based on its ve¬ 
locity, it can have no internal velocity structure. As a result, the 
kinetic energy below a certain mass scale ( actually the mass of an 
SPH particle and its neighbours) is not included in the calculation. 
We therefore stress that our simulation is lacking the kinetic energy 
that should be present at scales of less than ~ 2 OM 0 . The details 
of how individual stars form are thus not available from this simu¬ 
lation, and we must restrict ourselves to the large scale properties 
of star formation and the formation of clusters. 
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Figure 1. The panels show column density images from the simulations. The first panel shows the state of the gas after the crossing time (ter ~ 

and the following panel shows the gas at the free fall time. The remaining two images show how the cloud evolves to form a system of clusters. The maximum 

column density plotted is 0.12, 1.21, 2.90 and 4.86 g cm“^ respectively and the minimum column density is 10“^ g cm~^ in all the images. 


3 GENERAL EVOLUTION 

FigureQshows column density images from different points during 
the simulation and allows us to see clearly the evolution of the gas 
and regions of star formation. We see from the figure that the struc¬ 
ture of the cloud changes remarkably quickly. It starts as a churning 
network of gaseous filaments and within lOMyr (when the simula¬ 
tion was terminated) evolves into an ensemble of distinct clusters, 
by which time the gas has lost much of its early character. The fact 
that an unbound GMC can form stars and star clusters reinforces 
the predictions made in iciark & Bonnell| i2004) . 

The point at which the first bound objects condense out of 
the unbound flows occurs at roughly 2.4 Myr. This is roughly 
half the crossing time for the region (although some authors use 
ter = R/V instead of ter = ‘IRjV as is used he re). This time is 
consistent with the kinetic energy dissipation rate ilVIac Low et alJ 


ll998HStone eTailllQQSh and the formation of a turbulently domi¬ 
nated density structure JPadoan et all200lh . 

Rather than simply discuss the individual protoclusters that 
form (the sink particles) it makes sense here to discuss the bound 
groups of these protoclusters as well, which we will simple refer 
to here as ‘star formation centres’ or SFCs. The formation of the 
SFCs actually occurs very rapidly. The mass of the 16 most massive 
of these centres is shown as a function of time in figure]^ Within 
5Myr (or ^ 2.5 Myr after the onset of star formation) the 15 most 
massive SFCs all have masses greater than about 100 Mq, and are 
beginning to get to a size where there is good possibility of them 
forming massive stars (this will be discussed in section|3- 

A desirable feature of an initially unbound GMC is that cloud 
dispersal and star formation are occurring simultaneously. This re¬ 
moves the necessity for feedback mechanisms to disperse the cloud, 
or at the very least makes their task much easier. The timescale for 
star formation is thus comparable to the timescale for the cloud’s 
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Figure 2. Shown in the left hand panel are the masses of the 16 most massive SFCs and how they evolve in time. The horizontal long-dashed line gives the 
mass at which a SFC can form a star of > IOMq and the short dashed line shows the mass needed by a SFC to form a star of > 25M0 . This assumes the 
IMF given in section U] along with a star formation efficiency of 50% in the SFCs (see section Bl for a more detailed discussion). The solid line in the right 
hand panel shows the mass contained in the SFCs (gas particles -I- protocluster/sink particles). The dotted line in the plot shows the mass in stars greater than 
IOM 0 and the short-dashed line shows the mass in stars greater than 25M0 . The vertical long-dashed line denotes the point at which we estimate a supernova 
event to occur (note that tff= 4.7Myr). 


dispersal. The dynamics of an unbound cloud is therefore natu¬ 
rally in keeping with the recent observations that star formation 
and cloud dispersal occur in a few crossing times. There is also the 
added bonus that star formation efficiencies will be kept low, since 
most of the gas around a protocluster clump will be unbound to it 
and moving away. This prevents the material getting involved in the 
accretion once a SFC starts to form. 

Figure 0 shows the density distribution of the gas at three 
points in the simulation. The vertical dot-dashed line marks the 
original density of the cloud. Just before the first protocluster 
forms at 2.4 Myr we see that the most common density is roughly 
7 X 10“^^ gcm“® (the solid line curve), an order of magnitude 
higher than at the start of the simulation. Note however that very 
little material at this point is as dense as 7 x 10“^^ gcm“®, show¬ 
ing that the turbulence does n ot allow much material to get up to 
typical star forming densities iFalgarone, Phillips & Walkeill99ll : 

IPadoanI 1 99.51: IZinneckel^2002^. 

After 7 Myr the peak in the distribution falls back to roughly 
the starting density, however there is much more spread in the dis¬ 
tribution. This spread is controlled by two mechanisms. The high 
density tail increases as the SFCs grow by accretion and the sub¬ 
sequent rise in the potential energy. This causes yet more material 
to fall into the star forming regions. The low density tail increases 
since the cloud is freely expanding. By 13 Myr, only ~ Sfcr, we 
see that the majority of the gas has fallen to very low densities. 
By this point it is unlikely that observations of such a cloud would 
reveal much in the way of molecular gas and would instead only 
be visible as HI. The cloud can now be assumed to be ‘dispersed’. 
Even if the CMC fails to be a site of massive star formation, the dis¬ 
persal would still occur on a timescale consistent with Elmegreen’s 
(2000) observations. 

Note also from figure Q that the cloud contains cavities and 
dense regions of star formation. These are created in the simula¬ 
tion purely by the turbulence. This type of structure in star forming 


clouds is often attributed to the effects of high mass stellar feed¬ 
back, such as winds and supernovae, an d is thought to be the trig- 
ger for star formation in the region (e.g. lElmegreen & LadallQT^ 
Instead, we realise that turbulence can mimick these effects. Fur¬ 
thermore the cavities in the simulation woul d be easily ionise d by 
any high mass stars that form in the SFCs iPale et 'aHl2004 . We 
would then have a series of SFCs separated by a region of HII gas. 
Just as is found in the classic picture of triggered star formation. 


4 THE FORMATION OF STARS AND EXPECTED 
EFFICIENCY 

In this section we use some simple assumptions about the star for¬ 
mation that occurs in the SFCs to determine the numbers of high 
mass stars and the star formation efficiency that one might expect 
from the simulation. It is still beyond the capabilities of current 
computational resources to model the details of how individual 
stars form in a body of gas as large as a CMC. In the simula¬ 
tion presented here we cannot model any gas dynamics below the 
2 OM 0 scale. We can however give the reader a feel for the star for¬ 
mation that is present by using the results of previous simulations, 
along with some assumptions about the star formation efficiency 
and the form of the IMF. 

It has been shown from numerical simulations that star forma¬ 
tion occurs on roughly the local crossi ng time for the turbulence 
when the region is dynamically bound iBate et Ii] |2Q03l:lKk^ 
l200lh . On the small scales such as those represented by our proto¬ 
clusters, ~ O.lpc, the crossing time is of the order 10® years. We 
can therefore assume that all of our protoclusters form stars and that 
the star formation in our protoclusters takes place quickly, rapidly 
enough to be regarded here as instantaneous compared to the evo¬ 
lution of the whole GMC. 

The simulation presented also has no method of incorporating 
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logio ( Pgas(gcm-3) ) 

Figure 3. The plot shows the density distribution of the gas at three points 
in the simulation. The solid line follows the gas density just before the for¬ 
mation of the first protocluster, at 2.4 Myr (or ~ i ter). The dashed and 
dotted lines show the density distribution of the gas at a time of ~ 7 Myr 
and ~ 13.2 Myr respectively. The vertical dot dashed line shows the density 
of the gas at the beginning of the simulation. Note that the simulation only 
evolves a gas of pure molecular hydrogen. 


feedback into the GMC model. As is shown in figure|^in the right 
hand plot, the mass accreted into the SFCs gradually increases as 
the simulation progresses. At the point where the simulation is ter¬ 
minated, 30% of the GMC has been accreted by the protoclusters. It 
is unlikely that this value is representative of how much mass would 
actually he involved in the star formation hy this time, since feed¬ 
back mechanisms such as ionisation, winds and supemovae would 
seriously alter the amount of gas that would he available for ac¬ 
cretion into the SFCs. What is needed is an estimate of when one 
would expect the star formation process to he halted by feedback 
mechanisms. This requires some knowledge of the star formation 
taking place within the SFCs. 

We have already pointed out that the protoclusters in the sim¬ 
ulation group into large SFCs. From now on in the paper we will 
use the details of these regions, rather than the individual protoclus¬ 
ters, to assess the nature of the star formation in the GMC. TableQ 
gives the details of the SFCs after 9 Myr. The masses quoted for the 
SFCs in table Q includes all particles (SPH and protoclusters) that 
fall within the radius of the region. The gas particle component is 
however quite small, generally less than 20%. 

Although the star formation efficiency of GMCs is thought to 
be in the range of 1 to 10%, at the cluster level it is thought to be 
about 20 to 5 0% depending on the reg ion (for a discu ssion we point 
the reader to lLada & Ladall2003l and lKrou^l200 lli . In this paper 
we assume that the star formation efficiency in our SFCs is 50%, 
but will include a discussion about the case in which 100% of the 
mass is turned into stars. The assumed efficiency here is high but 
this is deliberate since it actually assumes as little as possible about 
the effect that the feedback mechanisms from the young stars are 
having on the accretion processes in the SFC. We will also assume 
that the IMF of the stellar population in the SFCs follow a two 
step power law form, dN oc wit h a = 1.5 for 0.08 < 

tti/Mq < 0.5 and a = 2.35 iSalDeteill955li for 0.5 < m/MQ < 


100. This IMF, in conjunction with our assumption that 50% of 
the mass of the SFCs is turned into stars, allows us to estimate the 
stellar population produced by the simulation. 

As already mentioned in the previous section, figure shows 
in the left hand plot how the mass of the 15 largest SFCs evolves 
with time. The horizontal lines mark the point at which high mass 
stars can form. From our IMF model, 15% of the mass should 
be contained in stars with masses greater that IOMq. Thus a 
IOMq star will be present provided that there is 10/0.15 =67Mq in 
the stellar population. Applying our assumed star formation effi¬ 
ciency of 50%, the SFCs must therefore have a mass of 134M0 if 
they are to harbour a IOMq star. The horizontal long-dashed line 
in the figure denotes the point at which the SFCs achieve this mass. 
Doing the same for 25M0 stars, which should comprise 7.7% of the 
stellar mass in our chosen IMF, we find that the SFCs need to con¬ 
tain 25/(0.077 X 0.5) = 650Mq if they are to contain a 25Mq star. 
This is represented by the horizontal short-dashed line in the figure. 

From our simple assumptions about the small scale efficiency 
and the form of the IMF, we can estimate at what point in the sim¬ 
ulation the star formation process will be disrupted by feedback 
mechanisms. From figure]^ we can estimate that the formation of 
IOMq stars would occur at about O.St/y (or at ^ 4 Myr). A star 
of mass 25Mq would form after ~ 1-lt// (or ~ 5 Myr). Since the 
mass of the SFCs is increasing fairly rapidly at this point, stars with 
even higher masses would be expected to be present shortly after 
this, within 0.5 Myr or so. It would thus appear that the GMC is 
able to get enough mass into the SFCs for them to be able to form 
a full stellar population within about 1 Myr. This is consistent with 
the observations of the small age spread in the stellar population of 
the Orion cluster iHillenbrand et all200ll) . 

Very rapidly after the first stars form we see that IOMq ob¬ 
jects will be present. This means that shortly after their formation, 
SFCs are going to contain ionising sources. Such stars are com¬ 
monly suggested to be responsible for controlling the star forma¬ 
tion efficiency by expelling the gas from the cluster in which they 
form (such as our SFCs), thus preventing the protostellar popula- 
tion from accretin g or preventing new stars from forming. However 
iDale et alJ j2004h have noted that the ionisation from these stars 
does not appear to significantly affect the accretion rate in the clus¬ 
ters. The clumpy/fractal nature of the gas at the centre of the cluster 
where the OB type stars are situated acts to shield vast regions of 
the cluster from ionisation. Rather than pushing through the dense 
material, the ionising photons just find the path with the least re¬ 
sistance out of the cluster. This is low density gas which would 
not normally be associated with protostellar accretion in the first 
place. Similarly, the gas structure may also prevent the winds from 
OB stars expelling gas from the cluster. It has been suggested that 
winds are able to escape via the fract al holes, witho ut imparting 
much momentum to the dense regions iHennin gll989ft . 

It is therefore not clear if ionisation or winds will be able to 
expel the gas from cluster, thus halting the star formation process. 
One mechanism that certainly will produce the desired effect is a 
supernova explosion. In fact it has been estimated that these events 
will not only remove the gas from a cluster, but also be able to 
disperse the natal GMC. Thus a high mass star’s death will defi¬ 
nitely mark the end of the star formation period in our cloud. Stars 
with masses greater than 25M0 have very short main-sequence 
lifetimes, of about 3-5 Myr, and we see from the figure that they 
form at about ~ 5 Myr after GMC formation. If we assume that a 
supernova event will occur at about ~ 4 Myr after the formation of 
the very high mass stars, then we estimate the first supernova event 
to occur at about 9 Myr or when the OB stars are ~ 4 Myr old. 
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Figure 4. Plotted are the positions of the SFCs large enough, by t = 9 Myr, to contain stars of mass greater than IOMq. The plot on the left shows the 
positions of the SFCs at t = 9Myr and the plot on the right shows how the SFCs will be positioned after 13Myr, assuming that a SN event expels the GMC gas 
and the SFCs continue on their original path. The circles show the size of the SFCs after the 4 Myr period, assuming they expand with their internal velocity 
dispersion. See section|^for a discussion. 


Assuming the supernova event will halt the star formation, 
we can now get an estimate of the star formation efficiency in 
the GMC. The vertical dashed line in figure denotes the point 
at which we might see the first SN event. At this time, ~ 0.1 to 
0.2 of the GMC’s mass is contained in the SFCs. However we have 
assumed up until now that the efficiency in the SFCs is not 100% 
hut 50%, therefore our estimate of the star formation efficiency in 
the GMC is roughly 5 to 10%. This is easily comparable to the 
expected efficiencies in GMCs by Elmegreen’s (2000) rapid cloud 
formation/dispersal model. 

The above analysis relied on a lot of assumptions about the 
nature of the star formation in the SFCs. In particular, it is guilty 
of invoking a star formation efficiency in the SFCs in order to de¬ 
termine the star formation efficiency of the cloud: one could argue 
that this is not entirely self-consistent. Here we redo the above anal¬ 
ysis but without the SFC efficiency assumption. If all the mass in 
the SFCs is used in forming a stellar population, then the mass re¬ 
quired by a SFC before a 25M0 star forms is 25/0.077 = 325M0 . 
The first SFC to achieve this mass does so at about 4 Myr, only 1 
Myr less than our previous estimate. Thus we can predict the su¬ 
pernova to occur at ~ 8 Myr. At this point in the simulation there 
is about 7 - 8% of the cloud incorporated in the SFCs. Thus the 
expected star formation efficiency is still less than 10%, suggesting 
that our analysis is not heavily dependent on the assumptions. 


5 DYNAMICAL EVOLUTION OF THE CLUSTERS AND 
THEIR RELATION TO OB ASSOCIATIONS 

As already discussed, the simulation produces a series of star for¬ 
mation centres (SFCs). Of these regions, 16 of them (see tabled 
are massive enough to contain a star of greater than 10 M 0 by the 
time at which we estimate a SN explosion will destroy the GMC. 
From figure Q we see that these SFCs are expanding as a group 
away from one another (in fact this is true of the GMC structure 
in general). Furthermore, the distance between SFCs is roughly 



distance from SFC centre of mass (pc) 

Figure 5. The figure shows the age of the protoclusters (note, not the SFCs- 
but the individual ‘sink particles’) as a function of distance. The ages are 
shown at a time of 13 Myr after the cloud formation. The distance shown is 
the y component from the centre of mass of all the protoclusters. The crosses 
denote the age based on when the protoclusters form and the hexagons 
points denotes the age based on when an individual protocluster can har¬ 
bour a star of greater than 10 M 0 . See sectionlslfor more details. It is clear 
from the distribution of points here that there is no discernable age spread 
with position in the region. The star forming centers are essentially coeval. 


lOpc after about 13Myr. Thus the group of clusters have the ap¬ 
pearance of an OB association, with the individual SFCs being OB 
subgroups that are expanding about some common point. In this 
section we examine the properties of the SFCs and compare them 
to the observations of OB associations. 
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Table 1. The table gives the properties of all the star formation cen¬ 
tres (SFCs) that would be expected to contain massive stars by t = 9 
Myr (see sectionl^for a discussion of this). AV is the internal SFC ve¬ 
locity dispersion and assumes the region is in virial equilibrium such 
that AV = {GMgpc/R)^^^- The crossing time is then calculated from 
ter = 2R/AV. Note that Mgrc is the total (gas -H stars) mass enclosed 
within R. 


SFC No. 

MsfC 

(M0) 

AV 

(km s“^) 

R 

pc 

ter 

Myr 

1 

1763 

3.27 

0.71 

0.42 

2 

761 

2.68 

0.45 

0.33 

3 

706 

3.24 

0.29 

0.17 

4 

624 

2.24 

0.53 

0.46 

5 

455 

2.37 

0.35 

0.29 

6 

362 

1.77 

0.50 

0.55 

7 

338 

1.77 

0.46 

0.51 

8 

329 

2.88 

0.17 

0.12 

9 

305 

1.68 

0.46 

0.54 

10 

212 

2.31 

0.17 

0.14 

11 

174 

1.37 

0.39 

0.56 

12 

151 

1.95 

0.17 

0.17 

13 

150 

1.95 

0.17 

0.17 

14 

144 

1.91 

0.17 

0.17 

15 

141 

1.63 

0.23 

0.27 

16 

140 

1.39 

0.31 

0.44 


In figure |4| the left hand panel shows the positions of the 
SFCs that are large enough to contain stars greater than lOM 0 at 
t = 9Myr, assuming a star formation efficiency of 50% and the IMF 
presented in section^ Their positions are plotted at the time we es¬ 
timate the SN event to start expelling gas from the cloud. We now 
assume that the SN event removes the gas from the GMC quickly 
enough such that the motions of objects have no time to adjust to 
the change in potential. We assume that this is true both at the scale 
of the SFC motions and at the smaller scale of the stellar motions 
inside the SFCs. The right hand panel shows the positions of the 
SFCs after a further 4 Myr, i.e. at t = 13 Myr, assuming that they 
continue on the path they had before the gas expulsion. The cir¬ 
cles denote the size of the SFCs at t = 13 Myr, which have been 
evaluated from r = rsN + AV x 4 Myr, where AV is the region’s 
internal velocity dispersion and rsN is the radius of the SFCs at the 
point when the SN event occurs (i.e. at ~ 9 Myr). This assumes 
that the star forming regions will disperse at roughly their internal 
velocity dispersion once the gas has been expelled. 

Figure|3clearly shows that the SFCs are expanding away from 
one another. By determining the average radius of the SFCs from 
their common centre of mass, both at t = 9 Myr and at t = 13 Myr, 
we can get an estimate of their expansion velocity from their dy¬ 
namical centre. At t = 9 Myr the average distance of the SFCs from 
their centre of mass is 18.4pc, while after 13 Myr it is 25pc. This 
corresponds to an expansion velocity (that is a 3-dimensional ve¬ 
locity) for the SFCs of 1.5kms“^. We compare this, for example, 
to the observed ex pansion of the OB subgroups in Per OB which 
is roughly 2kms~^ iFredrickll95(tlBlaauwll9^ . We also see that 
if the SFCs themselves are able to expand after the gas expulsion, 
they become an extended distribution of stars after only 13 Mys, as 
is shown by the circles in the figure. 

The stellar mass density is another important feature of OB 
associations. Generally, OB associations have a density of OB 
stars of roughly O.IMq pc~® (see the introduction for refer¬ 
ences). In this simulation at t = 9 Myr, the density of OB type 
stars is O.I 6 M 0 pc“®, and the mass density of stars > 25M0 is 


O.IM 0 pc“®. This is calculated by taking the total mass in the 
SFCs of stars of greater than the required mass type and dividing 
by the volume of region containing all the SFCs with these types of 
stars. After the system has had time to evolve for 4 Myr, the den¬ 
sities are 0.06 and O.O 4 M 0 pc“® for the OB type stars and those 
with masses > 25M0 respectively. Note these figures are based on 
the SFCs having a star formation efficiency of 50% and containing 
the IMF of stellar objects that was presented in section|4| 

We can compare this to the density of high mass stars in the 
SFCs at the point of the SN explosion. If we take the largest SFC , 
with mass 1763M0 and radius of 0.7pc, and assume again that 50% 
of this is contained in stars. Then the total mass in stars of mass 
greater than IOM 0 is 1763 x 0.15/2 = 132M0. The density of 
massive stars is then 132/0.7® = 384 M 0 pc“®. 

The turbulent flows are thus able to create a series of star form¬ 
ing regions that have roughly the same properties as those found in 
OB associations. Since the regions (the SFCs) are formed within 
large flows, the stars that form will have roughly the same mo- 
tion as the gas stream that formed them, potentially explaining why 
lBlaauwHl99lh finds that the gas surrounding OB subgroups is gen¬ 
erally moving with the group. 

Observations of some OB associations also indicate distinct 
age spreads between their subgroups. This generally takes the form 
of an age progression from one side of the association to the next. 
The ages are generally derived from where the very high mass stars 
turn off the main sequence, which is a much more reliable method 
than using pre-main-sequence (PMS) tracks of low mass objects. 
Also the high mass end is normally the only part of the mass spec¬ 
trum that is well established in OB associations. This age spread 
has been the motivation b ehind the triggered sequen tial star for¬ 
mation model developed bv lElmeg reMl&Mjtda^j^7^■ which is in 
turn motivated by the observations of ^Blaau^ll964^ . However, we 
note here that the Orion OB a ssociation exhibits no discernible age 
progression in the subgroups iBrown et alll999h . 

Does our simulation show a convincing age spread between 
the subgroups/SFCs? In figure we plot the age of the protoclus¬ 
ters (the sink particles that group together to form the SFCs) and 
their position from their common centre of mass. All the points 
are plotted at t = 13 Myr, with the positions being the y-direction 
in the simulation, since the SFCs are more spaced out in this di¬ 
rection. The ages are determined in two ways. The crosses denote 
the ages determined by when the protocluster first forms. The filled 
hexagons are determined by the time when the protoclusters reach 
a mass of 134M0 , the point at which a IOM 0 star can form. Note 
that in the previous sections we determined when IOM 0 stars could 
form based on the mass contained in an SFC , rather than its con¬ 
stituent protoclusters and gas. We are forced to use the individual 
protoclusters here since the SFCs are not coherent objects through¬ 
out the entire evolution of the simulation. 

We see clearly from figure that while a large range of ages 
exist at any particular distance from the centre of mass, no trend 
is present in the ages with distance. Thus our simulation predicts 
that the OB association would be essentially coeval. However this 
is just a symptom of our idealised initial conditions. The initial uni¬ 
form density sphere, with multiple Jeans masses, allows the entire 
cloud to proceed directly to star formation, via the dissipation of ki¬ 
netic energy. Since the turbulence is the same throughout the cloud, 
star formation occurs simultaneously in quite separate locations. If 
on the other hand our GMC needs to accum ulated in a large scale 
shock, as suggested bv IPringle et alJ i200ll) . then there would nat¬ 
urally be an age spread as the layer in which the GMC forms starts 
to grow. The most important point in this picture is that the whole 
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region would not be at the same density, but instead would have to 
evolve to star forming densities as the GMC accumulates. 


6 CONCLUSIONS 

The simulation presented in this paper highlights that GMCs need 
not be regarded as objects in virial equilibrium, or even bound, for 
them to be sites of star formation. Globally unbound GMCs can 
form stellar clusters very quickly, on roughly their crossing time. 
Furthermore, the unbound state of the cloud ensures that whole re¬ 
gion is also dispersing while it is forming stars. They are thus nat¬ 
urally transient features. This evolutionary picture of a cloud form¬ 
ing, producing a stellar population, and then dispersing has been 
shown by Elmegreen (2000) to be apparent in a number of inde¬ 
pendent observations. 

Using some simple assumptions about the form of the star for¬ 
mation in the star formation centres (SFCs) of our simulation, we 
have provided an estimate of the star formation efficiency in the 
GMC. At the point one would expect the first supernova events, 
we find that the star formation efficiency is about 5 - 10%. This 
assumes that the SFC environment in the simulation yields an effi¬ 
ciency of ~ 50%. Removing this assumption about the SFCs, and 
letting the SN event be the only control over the efficiency, we find 
that the cloud has a global star formation efficiency of ~ 7-8% (for 
our assumed IMF). 

We argue that unbound GMCs may provide a simple mech¬ 
anism f or forming OB a ssociations, a concept that was touched 
upon by I AmbartsumianI il95.5l [iQSSh . OB stars form at the cen¬ 
tre of a population of SFCs. These SFCs, which condense out of 
the unbound flows in the GMC are naturally expanding away from 
one another, as the positive energy disperses the cloud’s gas. Not 
only does the mechanism explain the OB association dynamics but 
it also explains the observed substructure, generally referred to as 
OB subgroups. Since the OB association is just a series of inde¬ 
pendently formed clusters, one would also expect the association 
to have the field star IMF. 
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